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To study a possible effect of transverse fluctuations of a stripe in a two-dimensional antiferromagnet 
on its charge dynamics, we identify elementary excitations of a weakly doped domain wall in the 
Hubbard model. Hartree-Fock numerics and analysis of fermion zero modes suggest that for U > 3t 
charged excitations are mobile holons, Q = 1, S — 0. Each holon resides on a kink in the position 
of the domain wall. We construct a simple model in which transverse stripe dynamics is induced 
solely by motion of the holons. In the absence of spin excitations (spinons, Q = 0, S = 1/2), stripe 
fluctuations do not suppress a tendency to form a global charge-density order. 



Physics of charge carriers confined to an antiphase do- 
main wall in a two-dimensional antiferromagnet (AF) is 
one of the most inipprtant and least—understood ingre- 
dients of the stripeBH-based approacha to high-TcpSiiper- 
conductivity. Both the lonewiistance interactionstffl and 
the short-distance physics§cl seem to be important for 
understanding even the basic phenomenology of stripes. 
GivenBQcl the existence of metallic stripes in an antifer- 
romagnet, their physics should not be far different from 
that of a more conventional one-dimensional electron gas. 
(1DEG). However, as suggested by model calculations ,t2l 
a metallic stripe with repulsive interactions placed in an 
AF environment may develop a spin gap. (In a conven- 
tional spin-rotation-invariant 1DEG a spin gap .can only 
appear in the presence of attractive interactions.!!!) With 
spin degrees of freedom gapped, coupling between adja- 
cent stripes could lead either to a global superconducting 
order- or to an insulating charge-density wave (CDW) 
state.0 In a 1DEG with rcpulsive_.interactions, CDW 
correlations diverge more strongly,EiJ and one would ex- 
pect them to dominate superconductivity. A way around 
this difficulty was suggestedEj by Kivelson, Fradkin and 
Emery (KFE), who noticed that transverse fluctuations 
of conducting stripes would induce a phase mismatch be- 
tween CDWs on neighboring stripes, effectively suppress- 
ing the interstripe coupling and the tendency to form a 
charge-ordered state. In contrast, such fluctuations do 
not suppress superconducting ordering, since the super- 
conducting phase is not spatially modulated. .— . 

Transverse dynamics of insulatingcj and metallictll 
stripes has been previously studied by Zaaiaen et al. Of 
a particular interest to us is the exampleo of a metal- 
lic stripe with a built-in coupling of charge motion and 
transverse fluctuations. The model describes an initially 
insulating stripe with a 4fcp CDW (1/2 electron per unit 
cell). Assuming that doped holes reside on rows adjacent 
to the site-centered stripe, Zaanen et al. found charge-1/2 
solitons, which shift the domain wall one lattice constant 
sideways (left/right), thus acquiring a transverse flavor. 

In this paper, we construct a similar model to study 
the effect of transverse fluctuations of a stripe on its 
charge dynamics. First, we employ Hartree-Fock (HF) 
numerics and an analysis of fermion zero modes to iden- 



tify elementary excitations of a lightly doped domain wall 
in the Hubbard model. Of course, such a texture is an 
excited, only locally stable configuration within the HF 
approximation. Nevertheless, it could be a good start- 
ing point for understanding the nature of globally stable 
metallic stripes. We find that a weakly-doped domain 
wall is &onrf-centered. At moderate and strong coupling, 
U > 3i, holes doped into such a wall create kinks in the 
transverse position of the stripe (wiggles) and form mo- 
bile solitons with spin S = and charge Q = 1 (holons, 
Fig. [I]). Unlike their ID counterparts, holons on a domain 
wall have a transverse degree of freedom, an isospin: a 
right (left) wiggle is ascribed isospin p = +l/2( — 1/2). 
this respect, they resemble the solitons of Zaanen et al 
Second, we construct an effective ID model of a stripe 
whose transverse fluctuations are associated solely with 
the motion of the holons, while spinons (kinks of a similar 
nature but with Q = 0, S = 1/2), are frozen out. In con- 
trast with the KFE arguments,E3 transverse fluctuations 
of a stripe in our model do not break phase coherence 
of CDWs on neighboring stripes. Therefore, they do not 
suppress the tendency to global CDW ordering, although 
the effect may be restored with the inclusion of spinons. 
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FIG. 1. Spin-charge separation on a bond-centered domain 
wall (dashed line) in an easy-axis AF (J z S> J x ,t). A hole 
decays into a spinon (Q — 0, S — 1/2) and a holon (Q = 1, 
S — 0) with isospins as shown (see text). Open symbols 
denote y-integrated spin and charge. 

Numerical mean-field solutions. We have solved self- 
consistently Hartree-Fock equations for a linearly polar- 
ized AF, (s x ) — (s y ) — 0, (sg) 7^ 0, for small and interme- 
diate interaction strengths U = 2 . . . 4t of the Hubbard 
model. At half-filling, we have found locally stable self- 
consistent solutions with a bond-centered domain wall. 
A site-centered wall has a slightly higher energy. A re- 
duction of staggered magnetization on the wall induces 
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a ID band in the middle of the Hubbard gap A (Fig. ||). 
Fermions confined to a site-centered wall have a Dirac 
spectrum E = ±v\k x — ir/2\ and acquire a small mass 
(gap) A < A in the case of a bond-centered wall. This 
gap arises because an electron moving along the wall feels 
a nonzero x-staggered magnetization (vanishing by sym- 
metry for a site-centered wall). At half- filling, the Fermi 
energy lies between the two midgap bands (Ep = U/2). 
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FIG. 2. Left: midgap one-particle spectrum E(k x ) of an 
AF with a straight domain wall. Self-consistent solution of 
Hartree-Fock equations. U = 2.5t, 48 x 27 sites. Dashed line: 
low-energy approximation, Eq. (El). States outside the gap 
are not shown for the site-centered stripe. Right: 4 Fermi 
patches in the Brillouin zone, Eq. (^|), shown with eigenvalues 
( T 3! r 3)- Dashed line: Fermi surface of the noninteracting 
system. 

A bond-centered domain wall with a wiggle is similar 
to an AF chain with a kink in x-staggered magnetization 
(Fig. ||). As in the model of polyacetyleneJi3 a pair of 
degenerate localized states (zero modes) appears in the 
middle of the smaller gap. A wall with a kink doped with 
a single hole has a S = 0, Q = 1 soliton at the wiggle 
(Fig. |j). In view of translational invariance, the charged 
defect can move along the wall. We thus identify it with 
a holon of the large- J z limit (Fig. |l|). At U = 2 . . .3t 
(rather weak coupling) , helons are cigar-shaped like spin 
bags of Schrieffer et al.^3 a consequence of the Fermi 
surface nesting. Orientation of a holon along one or the 
other lattice diagonal is correlated with the direction of 
the wiggle (the transverse flavor). We have verified that 
holons are the lowest-energy states of doped charges on 
an empty wall for U > 3i (at the HF level). At a weaker 
coupling, they bind ito bipolarons with Q = 2, S = 0, 
and zero isospin. 
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FIG. 3. Left: Bond-centered stripe with a wiggle as a su- 
perposition of a site-centered domain wall (black arrows) and 
a ID AF chain with a kink (open arrows). Right: Holons with 
alternating isospins form a cite-centered stripe. 
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FIG. 4. Staggered magnetization m(r) = ( — l) x+y (s z (r)} 
and hole density 1 — n(r) in a Hartree-Fock calculation at 
U — 3i, 24 x 24 sites. A bond-centered wall with 2 wiggles 
and 2 doped holes. 

Continuum Hartree-Fock approximation. HF equa- 
tions for a linearly polarized Hubbard AF read 

- t J2 i>~s{v') + U[dn(r) - s m(r)ty,(r) = E^ s (r), 

<rr'} 

(1) 

where m(r) = ( — l) x+y (s z (r)) is the staggered magneti- 
zation and s is the staggered spin index. In what follows, 
density fluctuations Sn(r) will be neglectedEZl to restore 
charge-conjugation symmetry. Eq. (PD) can be cast in a 
matrix form using Pauli matrices {<ri}. Staggered spin s 
is an eigenvalue of 03, the hopping term is proportional 
to <7i . Most importantly, zero modes of Eq. (Q) are in- 
variant under charge conjugation tp(r) — > <72^*(r). If 
such modes are present, the system contains solitons-with 
Q = 0, S = 1/2 (spinons) or Q= 1,5 = (holons). E3 

We have determined (see, e.g., Fig. g, left) that low- 
energy midgap states induced by an empty domain wall 
live at one of the four "Fermi patches" with lattice mo- 



menta Ifor 



tt/2 (Fig. 0, right). Introducing 



smoothly varying amplitudes of a fermion wavefunction, 
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(2) 



adds two more indices, a = sgnfc^ and /i = sgnfc y . Ac- 
cordingly, we preserve only those Fourier components of 
magnetization which connect the Fermi patches: 



; (r)) ~Xia=0 Eu= 



(r)e 



With the new indices come two more mutually commut- 
ing sets of Pauli matrices, {rf} and {r*}, i — 1,2,3. In 
terms of these, k « f(rf,r|), (-If = rf, (-1)* = rf. 
The HF Hamiltonian becomes 

Hhf — —1itao\T^d x — 2itao\T^d y — Ua 3 m(r), (3) 
to = mn + JnoiTf + m 10 T% + m 00 T^Tf . 

Only TOn(r), staggered magnetization proper, survives 
in the bulk inducing the Hubbard gap A = J7|toh(oo)|. 

On a straight domain wall in the x direction, tooi = 
tooo = 0, while mirj 7^ 0. The spectrum of midgap states 
depends on the symmetry of the wall. If the stripe is 
site-centered (bond-centered), wall fermions have a gap- 
less (gapped) spectrum. The absence of a gap can be 
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ascertained by finding a zero mode at the momentum 

p x = -id x = 0: 



dip U 
a2 Ty = 2Ta [T * mil{y) 



■iT$m 10 (y)]il)(y). 



(4) 



To zeroth order in mio, solutions of Eq. (Q) are eigen- 
states of (72, Tg and, e.g., rf (each zero mode comes 
from a single Fermi patch) , giving .a. total of 8 linearly 
independent zero modes. As usual,E3 only half of these 
solutions [those with o~%t\ mn(+oo) < 0] are localized on 
the wall, so that there are 4 zero modes. Remarkably, in 
addition to the usual twofold spin degeneracy, there is 
another spin-like degree of freedom, which will prove to 
be the transverse flavor. The origin of isospin (at weak 
coupling) is thus exposed: compared to a ID chain, there 
are twice as many "Fermi points" on a straight domain 
wall in 2D — see Fig. g, right. 

The difference between one-particle spectra of site and 
bond-centered walls arises in the first order in mio. 
Eq. (^) has four zero modes if ra^p is an odd function 
of y, i.e., for a site-centered wall.Ej On a bond-centered 
wall, mio(y) is even and fermions have a gap (Fig. |J): 

E(0) = ±U{m x ) =±U [ m 10 (y)^{yW(y)dy. 



The one-particle spectrum E(p x ) on a straight site- 
centered stripe can be determined approximately by 
starting with p x = [Eq. ([I])] and treating the first term 
in Eq. (|J) perturbatively. In the limit p x — > 0, states 
outside the main gap can be neglected, which reduces 
the Hilbert space to the four zero modes [Eq. (|J)]. By 
using the degenerate perturbation theory, we find a Dirac 
spectrum (dashed lines in Fig. 0, left): 



E(p x ) ~ ±vp x , 



2ta(rl 



(5) 



Domain-wall holons at U <C t. As illustrated in Fig. |^, 
magnetization on a frond-centered wall with a wiggle 
can be obtained by superimposing m(r) of a straight 
s^e-centered wall and that of a spin chain with a kink 
in a;-staggered magnetization. Away from the wiggle, 
moo(r) = TOoi(r) = 0. To simplify the discussion, we will 
neglect these components altogether. Decompose m(r) 
into an x-independent part and the rest: 

m(r) = m (0) (y) + m (1) (r), 
m^°\—y) — —m^(y), mP^ (±oo, — y) — mP' (±oo, y). 

The Hamiltonian (||) can now be split in two parts: 

= ~2ita a.rldy - Ua 3 [m£ } (y) + (y)rf ] , (6) 
H™. = ~2ita ayrtd x - Ua 3 [mft (r) + m$ (r)rf ] . (7) 

As shown above, the "transverse part" (||) has 4 zero 

modes for each p x . Within this Hilbert space, fl^F 
scribes right and left-moving fermions with spin, which 
see a staggered magnetization 

(mi (a;)) = J dy v){y)[m w (r) + mn(r)rf ]u(y), 



where u{y) is a zero mode (§ of Eq. (§). The 
midgap fermion band acquires a gap of its own, A = 
[/|(mi(oo))| < A, with two zero modes (one for each 
spin) inside this smaller gap. "Longitudinal" wavefunc- 
tions of the two zero modes satisfy the equation 
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dip(x) 
dx 



U 
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r£(mi(x))Tp(x). 



(8) 



The existence of two holon flavors can now be deduced 
from Eqns. (Q) and (S). The zero modes have a finite 
norm only if 

<T2T$m$(+O0) < 0, <T27f <mi(+oo))/(7f) < 0. 
It follows then that the product of eigenvalues 



T x v _ 



sgn[m (0) (y = +oo) (mi (a; = +oo)) (rf)] (9) 

can be identified with the holon isospin 2p. This can be 
seen by extrapolating Eq. (|]) to larger values of U, which 
reduces the size of holons. We have (rf ) = ((— l) v ) = 
(— l) yo , where yo is the row number of the chain in Fig. [|. 
According to Eq. (g), if rf r| = +1, spins on the chain 
and to the right (left) of the wiggle are an extension of 
the upper (lower) AF domain, as for the p = +1/2 wig- 
gle in Fig. Thus, p — rfr|/2. This identification is 
consistent with numerical HF solutions (Fig. ||), where 
T 3 T 3 = s g n ^r sgnfcj, can be inferred from the orientation 
of a holon. 

Effective holon Hamiltonian. Large- J z cartoons (e.g., 
Figs. [I] and [sj) suggest that any configuration of a do- 
main wall is uniquely represented by an interface drawn 
through integer (if a hole is present) or half-integer (bond, 
no hole) lattice points. In the absence of overhangs, the 
imaginary-time interface dynamics due to spin exchanges 
and hole hops can be described by an SOS-type model. 
(Locality in time assumes-that coupling to the bulk spin 
waves can be neglected.EJ) In the dilute limit, mobile 
single-particle excitations are spinons and holons, each 
equipped with an isospin 1/2. When spinons are frozen 
out, EH the only remaining excitations are holons, and the 
SOS model reduces to that of (iso)spin-l/2 fermions with 
the short-range Hamiltonian 



H 



E[(- 



<J X s, 



h.C. 



J Z Si &i 



(10) 



+U n iT mi +Vm n l+1 -r , Jz ^ a i+1J 

Here m = n i] + n i[ l& the holon-number operator, 
m p = tpjpipip, and isospin operators are sf = (n^—nii)/2 
and sf = iptftpii- Taking U — * +oo excludes double oc- 
cupancies, while isospin exchange J x allows two holons 
on one side of the stripe to hop to the other side. By con- 
struction, we do not expect isospin-rotation invariance; 
in general, J z ^ J x . 

Holons can be transported between different stripes 
without disturbing the AF order only in pairs (with op- 
posite isospiiis). In contrast to the model of a quarter- 
filled chairO with the same Hamiltonian (10), the field 
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tp^ describes Q = 1 particles, and the appropriate pairing 
operator V'l+ii nas charge Q = 2. This operator is a 
combination of the (iso)singlet (SP) and an (iso)triplet 
(TPq) pairing operators, 



O sp , O- 



TPn 



and an instability in either of these channels could lead 
to global (Cooper-pair) superconductivity, even though 
isospin is local to a given stripe. 

To analyze this scenario ,_aie_ have performed a stan- 
dard weak-coupling analysisEilEj at a generic density. Be- 
cause a stripe horizontal on average is invariant under 
isospin reflection, isospin and charge degrees of freedom 
separators The scaling in each sector is determined by 
the usual constants K s , K c ; for repulsive interactions 
K s > K^ 1 > 1. In the absence of the isospin gap, both 
SP and T-Eo correlation functions have the temperature 
exponent^! fisp = Mtp = K~ l + K s — 2 > 0, i.e., there 
is no divergence in either channel. Unlike in a more con- 
ventional IDEGjij where the spin-rotation invariance re- 
quires J z — J x and another mechanism {e.g., the "spin- 
gap proximity effect"E3) is needed to develop the spin 
gap, here an isospin gap arises naturally for J z > \J X \ 
(an easy-axis anisotropy). At low temperatures, such a 
system is in the Luther-Emery phase, the TPo component 
freezes out, while the isosinglet pairing exponent becomes 
£tg P = K~ x — 2; it diverges for K c > 1/2. Isospins are 
ordered at T = and the stripe becomes, on average, 
site-centered (Fig. ||, right). 

As usual, the more divergent (/xcdw = K c — 2 < /xsp) 
CDW correlations compete with the superconducting or- 
der. In oux_ holon-only model (Jl0|) , transverse stripe 
fluctuations^} never suppress this instability. CDWs on 
neighboring stripes remain phase-coherent not only in 
the isospin-ordered phase (where stripe fluctuations are 
suppressed by the isospin gap), but also in the presence 
of fully developed transverse fluctuations. 

Large- J z cartoons indicate that the length of a stripe 
seen by a hole increases in the presence of spinons (cf. 
Ref. [l2]). This is why we believe that in the full inter- 
face model (with both spinons and holons present, as in 
Fig. |l|) the transverse fluctuations may reduce the coher- 
ence of CDWs on neighboring stripes. However, other 
properties of this model require a separate study. 
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